function [output1,output2] = comeds2D(input)
% function [output1,output2] = comeds2D(input)
% This function is the driving function for the three different
% coars mesh solvers implemented
% The input contains everything that each solver needs
%   ...maybe a data structure??


% POWER ITERATION
% outputPI = solverPI(input)

% NEWTON
% outputNT  = solverNT(input)

% NEWTON-GMRES
% outputGM  = solverNG(input)


% return output